Quantum optical effective- medium theory for loss-compensated metamaterials 
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A central aim in metamaterial research is to engineer sub-wavelength unit cells that give rise to 
desired effect ive-medium properties and parameters, such as a negative refractive index. Ideally one 
can disregard the details of the unit cell and employ the effective description instead. A popular 
strategy to compensate for the inevitable losses in metallic components of metamaterials is to add 
optical gain material. Here we study the quantum optics of such loss-compensated metamaterials 
at frequencies for which effective parameters can be unambiguously determined. We demonstrate 
that the usual effective parameters are insufficient to describe the propagation of quantum states of 
light. Furthermore, we propose a quantum-optical effective- medium theory instead and show that 
it correctly predicts the properties of the light emerging from loss-compensated metamaterials. 
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Noble metals are an important ingredient of metamateri- 
als, i.e. man-made structures composed of strongly sub- 
wavelength unit cells, with effective dielectric parameters 
often not occurring in nature, such as a negative refrac- 
tive index [l], Q . Metals are also inherently lossy, and for 
many applications it is naturally desirable to have less 
loss. Complementary strategies are to replace the met- 
als by other material [3] or to compensate for the metal 
loss 4-8]. As an important branch of active plasmon- 
ics |9j , active loss compensation with the use of gain ma- 
terial has already proved experimentally successful, for 
example in surface plasmon polariton propagation [Tq| 
and in metamaterials [Tlj . 

The constituents and geometry of a unit cell can be 
complicated and interesting, but they are designed to 
allow the effective description of the metamaterial as 
a homogeneous medium. For subwavelength unit cells, 
unique effective dielectric parameters can be identified, 
independent of the method used to retrieve them. Then 
the common belief is that experiments at the operating 
frequency do not reveal information about the unit cell 
beyond the usual effective refractive index. For larger 
unit cells, the extraction of effective material parameters 
becomes less meaningful, not just depending on the unit 
cell, but also for example on the size of the metamaterial. 



In recent years the quantum optics of attenuating 12 



17j and amplifying [18R22| dielectric media was devel- 
oped, where optical modes are described as open quan- 
tum systems. There are important similarities with clas- 
sical optics, for example the classical Green function 
plays a central role also in quantum optics, but new 
is that both with loss and gain there is quantum noise 
associated. This leads to the important question, not 
addressed before to our knowledge, whether the usual 
effective-medium theories for metamaterials that work 
well in classical optics can also be used in quantum op- 
tics. 
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FIG. 1: (Color) Sketch of the setup: Quantum states \L,R) of 
normally incident identically linearly polarized light probe the 
metamaterial from both sides, and the output on the right is 
analyzed with balanced homodyne detection. The multilayer 
metamaterial in air has alternating layers with thicknesses 
<i a ,b and volume fractions p a ,b = d a ,b/d, with d = d a +dh- The 
two outermost layers have widths d a /2, making the metama- 
terial finite periodic with M left-right symmetric unit cells. 
In Figs. [TH3] we choose d = 0.2c/cjo and M = 5. The suscep- 
tibilities s a ,b((^) as in Eq. (pQ) describe either loss or gain. 



In this Letter we consider the propagation of quantum 
states of light through simple types of metamaterials, for 
which well-established methods to retrieve effective pa- 
rameters agree very well. We show that in particular for 
loss-compensated metamaterials the output states can be 
considerably different than for their homogenized coun- 
terparts, when employing the usual effective- medium the- 
ories. In response to this observation, we propose an ac- 
curate quantum optical effective-medium theory. 

Loss- compensated multilayer metamaterial. — We con- 
sider loss-loss, gain-gain, and loss-gain periodic multi- 
layers, the latter being loss-compensated metamaterials, 
see Fig. [lj Linear gain occurs when pumping the gain 
medium below the lasing threshold [23]. For definiteness 
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we describe the dielectric response of the layers by single- 
resonance Lorentz oscillator models [l9[ 



e(u>) = 1 + 



/ N down 1 



up 



(i) 



which describe media of two-level atoms near a resonance 
frequency co>o, with N up < Adown and e\ = Im(e) > 
for lossy media, while gain media have inverted popula- 
tions iV U p > A/down and hence s\ < 0. Eq. (pp) satisfies 
Kramer s-Kronig relations, and causality is essential for a 
consistent quantum optical description. One can param- 
eterize the populations occurring in Eq. (pQ) by a thermal 
distribution iVth(cj, \T\) = [ex.-p(huj/kB\T\) — with an 
effective temperature T that equals N up / (Adown — ^up) > 
for lossy and N dovm /(N up - Adown) < for amplifying 
media [19| ; T = K corresponds to no excitations in the 
former and complete inversion in the latter. 

Effective parameters, two methods. — We consider non- 
magnetic layers, giving a nonmagnetic effective medium 
with /i e ff = 1. To determine the effective dielectric func- 
tion, we compare two well-established methods: first, we 
take the volume average £ e ff,i = Pa^a + Pb^b, the ar- 
gument being that the optical wavelength is so much 
larger than a unit cell that it does not probe its in- 
ternal structure, but only its average response. As our 
second method, we determine the effective permittivity 
£eff,2 from the complex amplitude reflection and trans- 
mission coefficients r, t of the metamaterial, as proposed 
by Smith et al. [lif. We determined £ e ff,i and £ e ff,2 for 
all structures considered in this work, and as expected 
since d <C Ao for all of them, we find that they agree very 
well, as would other homogenization procedures [25l427j. 
Thus we denote "the" effective dielectric function by s e ff- 

Challenge in quantum optics. — For the propagation 
of quantum states of light through the metamaterial, we 
need to take into account that both with loss 12l4l7| 
and with gain TH-IHj there is associated quantum noise. 
We do this by generalizing the formalism by Gruner and 
Welsch [13| to structures with not only loss but also gain. 
In this input-output formalism, the output operators are 
expressed as linear combinations of the input operators 
plus quantum noise terms Fl,r- 



&L,out 
&R,out 



A 



(2) 



Here the 2x2 matrix A is the same that in classical optics 
connects input and output field amplitudes. The Fl,r(^) 
on the other hand have no classical analogues, and repre- 
sent linear combinations of the quantum noise operators 
associated with all layers (see below). The challenge for 
effective-medium theories in quantum optics is to predict 
observable effects of this quantum noise well. 

Homodyne detection of output light. — Let us now con- 
sider an experiment in which we probe our metamaterial 
with quantum states of light, and predict the output as 



measured with a balanced homodyne detector [23, l29j . 
see Fig. [TJ As input states, we take continuous-wave 
squeezed vacuum states of light 0. In quantum op- 
tics we need to take into account that the multilayer has 
two input ports, and to start with good squeezing, we 
prepare squeezed vacuum states in both. The squeezed 
state |L) = <Sl|0) (|R) = <Sr|0)) characterizes modes im- 
pinging on the left (right), with the squeezing operator 



<S L exp 



jy dcj [f£ Q>L,m(w) «L,in(2^ - Lj) - H.C.] j . 

(3) 

Here, £l = |£L|e 1( ^ L is the squeezing parameter and <2L,in 
are annihilation operators, see Fig. [TJ The analogous 
formula for <Sr features aR 5 i n (o;) and £r = |^R|e 1( ^ R . 

The homodyne detection involves a local-oscillator 
field with frequency cjlo and phase 0lo- The differ- 
ence photocount at both detectors is proportional to the 
output quadrature variance [3, l2lj 
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2Re[<oJ ut ,aL>e 2i * LO ] 



(4) 



where we wrote a ou t for aR )OU t(^Lo), and use short-hand 
notation (C, D) for (CD) - (C)(D), with (. . .) the quan- 
tum expectation value. In the time domain, the operator 
E(^LO,t) is given by a RiOUt (t) exp[i(u LO t- </>lo -tt/2)] + 
H.c This quadrature component of the light is called 
squeezed if the variance ([AE^lo? ^lo)] 2 ) is less than 
unity [1H, the value for the vacuum state in free space. 

Usual effective-medium theory in quantum optics. — In 
the usual effective-medium theory, the optical properties 
of a nonmagnetic metamaterial can be completely de- 
scribed in terms of the complex- valued effective refractive 
index n e ff = V^eff = ^eff + i^eff and the overall geome- 
try (here: thickness). The usual effective- medium theory 
when applied in quantum optics amounts to replacing in 
Eq. ([2]) the matrix A for the multilayer by another 2x2 
matrix AeS that describes transmission and reflection of 
light by a single slab with thickness Md and complex ef- 
fective refractive index h e n = n e ^ + ift e ff- Thus we have 
An exp(iujMd/c)A22 ^eff, and A12 = A21 t e ff, 
with t e ff and r e ff the usual transmission and reflection 
amplitudes for such a slab. Correspondingly, we replace 
the quantum noise operators El,r of the multilayer by 
the known operators EL,R,eff for a single slab with refrac- 
tive index h e n [l3|]. Important for the homodyne signal 
is the zero-average noise operator F G s = ER, e ff(^), that 
depends on the effective noise-field operators (p e s(z,Lj) 
everywhere in the homogenized slab, 



"eff 



|g e ff,l|^ ^eff 

4c n e ff 



dzf(z)(p eS (z). (5) 



Here, we suppressed the frequency dependence and wrote 
/ = /+ + /_ with f±(z,cj) = (n e ff ± 1) exp(±in eff cjz/c). 
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FIG. 2: (Color) Output quadrature variance of Eq. (|4]) as 
measured with balanced homodyne detection, for squeezed 
vacuum input states with |£l| = |£r| = 0.2 and 0l = 20lo — 
5 and 0r = 20lo — 2. Metamaterial dielectric parameters 
in Eq. (PQ): cj 0a ,b = u>o, u) pa /u = 0.3, j a /u)o = 0.1, and 
<^pb/<^o = 0.25 and jb/^o = 0.15. Effective temperatures 
T a = Tb = 0. Panel (a): loss-loss multilayer. Panel (b): gain- 
gain multilayer. Red dashed lines (usual effective medium 
theory) overlap blue solid lines (exact multilayer theory). 

The commutator [<£ e ff(z, gj), (p\^{z' ^ u')} equals 
sgn[e e f[ i i(uj)]S(z — z')S(uj — u/), and thus depends 
on the effective medium being lossy or amplifying. 

The electric-field quadrature variance Eq. ((4]) as 
measured with homodyne detection will depend on 
the quantum noise through the expectation value of 
(Fq S (a;).Feff (^O)j which f° r thermal noise becomes 

{JVthfl[e«ff ) i]-(JVth+l)e[-e«ff ) i]}(l-|r e ff| 2 -|teffr)<y(w-w / )- 

(6) 

Since (1 — |r e ff| 2 — |^eff| 2 ) is positive for net-loss and neg- 
ative for net-gain media, this expectation value ((6]) is al- 
ways > 0. In particular, the usual effective medium the- 
ory predicts that quantum noise can be neglected at ex- 
act loss compensation, £ e ff,i = 0. Notice also that for the 
limit T = OK, Eq. (O predicts (F^(u)F eS (u')) to van- 
ish for effectively lossy media, whereas the corresponding 
value for net gain media stays finite. 

Loss-loss and gain-gain multilayers. — As a first test of 
the usual effective- medium theory in quantum optics, we 
present in Fig. [2] the output quadrature variance for loss- 
loss and gain-gain metamaterials, all at T — 0. As the 
figure shows, the usual effective medium theory predicts 
the homodyne signal well, for both multilayers. For the 
loss-loss multilayer at zero temperature quantum noise 
can be neglected so the agreement only confirms that 
A can be replaced by v4 e ff, as we know from classical 
optics. By contrast, for gain-gain multilayers, the effects 
of quantum noise in the output cannot be neglected at 
zero temperature, but the figure illustrates that here also 
the usual effective- medium theory is accurate. The same 
holds true (not shown) for the structures of Fig. [2] at 
finite temperatures T a = Tb ^ 0. 

However, if due to pumping the effective temperatures 
T a and Tb are different in both layer types, then the usual 
effective- medium theory seems incomplete, since then it 
is not clear what effective distribution N g s(<jj) to use in 
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FIG. 3: (Color) Loss-compensated metamaterial with geom- 
etry of Fig. [U with lossy a-layers and amplifying b-layers. 
Material parameters and input states as in Fig. [2] Panel (a): 
frequency dependence of £ e ff,i, showing net loss around cjo 
and net gain (yellow) elsewhere. Exact loss compensation 
occurs at 0.77ujo and 1.30cjo- (b): Output quadrature vari- 
ances predicted with the usual effective- medium theory (red 
dashed line), and with the quantum optical effective medium 
(QOEM) theory (dash-dotted green line); the latter overlaps 
the exact multilayer calculation (solid blue line), (c): corre- 
sponding effective noise photon distribution N e f£ of Eq. (0). 

Eq. (|6]). One would need to postulate it, and if indeed e e ^ 
is to give a complete effective description of the medium, 
then one has little choice but to postulate that N e ^ equals 
the volume average 5^j= a bft^ th ( Co '' l^il)- However, this 
does not give accurate predictions, as we discuss below. 

Breakdown of usual effective- medium theory. — For 
loss-compensated metamaterials, we see in Fig. [3f b) 
already for T a = Tb = a clear failure of the usual 
effective-medium theory, which predicts the output light 
to be squeezed for almost all values of cjlo near the mate- 
rial resonance co>o, whereas the full multilayer calculation 
shows that actually the output state will not at all be 
squeezed in the whole range of c^lo/^o between 0.76 and 
1.21. The output of a loss-compensated material is nois- 
ier than of a homogeneous slab with the same n e ff. 

Quantum- optical effective medium theory. — We will 
now derive a quantum-optical effective medium (QOEM) 
theory that gives accurate predictions also for loss- 
compensated metamaterials. Our approach is to distill 
solely from a unit cell not only the usual £ e ff(^), but also 
an effective noise photon distribution iV e ff(u;). 

Applying the general multilayer theory to one 
unit cell only, we compute the noise contribution 
(^unit(^)^unit(k/)) t° the output variance. The use of 
the general theory ensures that we do not introduce the 
unphysical noise compensation that we found above to 
be inherent in the usual effective-medium description. 
Since the unit cell is much smaller than an optical wave- 
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length, it suffices to Taylor expand the average to first 
order in the layer thicknesses d a ^. Analogous to Eq. (j6j), 
we equate (Fj nit (cj)F unit (a/)) with {iV e ff0[£ e ff,i] - (N e s + 

uct of the yet unknown iV e ff and the net gain or loss of 
a unit cell, the latter computed in the usual effective- 
medium theory and also to first order in the thickness. 
Solving for iV e ff, we obtain as a main result the effective 
noise photon distribution 

iV eff H = -i + 1 Cj[2N th (u, \Tj\) + 1], (7) 

j=a,b 

in terms of the dielectric parameters of the unit cell 



Thus, our QOEM theory predicts that the quantum noise 
contribution (F^(u)F e ^(uj f )) to the output variance of a 
multilayer metamaterial consisting of arbitrarily many 
unit cells is given by Eq. (J6j) , but with the thermal distri- 
bution iVth replaced by the effective distribution N e s of 
Eq. (0. Fig.[3jb) illustrates the very good agreement be- 
tween the output variances computed with the exact the- 
ory and with the QOEM theory for a loss-compensated 
metamaterial. N G s grows when the loss is more exactly 
compensated by gain [smaller |£ e ff,i|> see Fig. G^a)] or 
when the same value |£ e ff,i| results from compensating 
more loss with more gain (larger When loss is ex- 

actly compensated N e ^ even diverges, while the output 
variance stays finite, compare Figs. [3ja,b,c). 

The amount by which the sum J2j =a b Cj exceeds unity 
is a measure of how much loss is compensated by how 
much gain, which is information beyond £ e ff (^0- For loss- 
loss and gain-gain metamaterials the ( a and £b add up 
to unity. Only in case T a = Tb, as for a passive lossy 
metamaterial at room temperature, we then find that 
iVeff reduces to the thermal distribution Nth- Thus our 
QOEM theory explains our earlier observation that the 
usual effective-medium theory worked well in Fig. O 

Conclusions and outlook. — We have shown that the 
usual effective-medium theory breaks down in quantum 
optics, and that the effective dielectric function e G Q does 
not always suffice to describe the output even of perfectly 
homogenizable metamaterials. In particular for loss- 
compensated metamaterials, we showed that balanced 
homodyne detection can be used to 'open the effective- 
index blackbox' and to detect in which of two metamate- 
rials with the same n e ff most loss compensation occurs, 
because loss compensation is not conveniently accompa- 
nied by quantum noise compensation. We showed that 
on this point the usual effective-medium theory fails in 
quantum optics, and also fails without loss compensa- 
tion when different effective temperatures coexist due to 
pumping. 



Moreover, we proposed and tested a quantum-optical 
effective medium (QOEM) theory that works for loss- 
compensated and for gain-gain and loss-loss metamateri- 
als alike. The effective dielectric function is the usual one, 
but it is complemented by an effective noise distribution 
N e s(uj) which involves dielectric properties beyond £ e ff, 
namely the positive loss and gain parameters Ca,b- These 
central results of Eqs. (|7|8j) can be generalized without 
difficulty to more than two layers per unit cell or contin- 
uously varying e(z). The formulae do not depend explic- 
itly on the unit cell being layered, which suggests their 
validity for a much larger class of metamaterials. 

Our findings are not just relevant for homodyne detec- 
tion of squeezed light, but for the optics of metamateri- 
als in general, for example when employing superlenses 
based on loss-compensated multilayers as proposed in 
Ref. y to focus single- or few-photon pulses of light. The 
quantum noise grows with thickness of the metamaterial 
and is independent of the optical input states, and thus is 
relatively more important for weaker input fields (fewer 
photons) and negligible for bright light. But whatever in- 
put state is used, there will be additional output variance 
due to the quantum noise in pumped metamaterials, un- 
derestimated by the usual effective-medium theory, but 
well predicted by our QOEM theory. Analogous to meta- 
materials, it would be interesting to study the quantum 
optics of loss-compensated slow-light media [HI . 

In summary, we identified a fundamental relation be- 
tween loss compensation and quantum noise, constitut- 
ing a breakdown of standard effective-medium theory, 
that should be taken into account for few-photon meta- 
material applications. We expect our quantum optical 
effective-medium theory to become a useful tool for the 
design and use of active metamaterials. 
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